Simple trace criterion for classification of multilayers 



L. L. Sanchez-Soto, J. J. Monzon, T. Yonte 
Departamento de Optica, Facultad de Ciencias Fisicas, 
Universidad Complutense, 28040 Madrid, Spaii^ 



o 
o 

(N 
< 



> 
00 
(N 
O 
00 

O 



Oh: 



X 



J. F. Carinena 
Departamento de Fisica Teorica, Facultad de Ciencias, 
Umversidad de Zaragoza, 50009 Zaragoza, Spam 
(Dated: February 1, 2008) 

The action of any lossless multilayer is described by a transfer matrix that can be factorized in 
terms of three basic matrices. We introduce a simple trace criterion that classifies multilayers in 
three classes with properties closely related with one (and only one) of these three basic matrices. 

PACS numbers: 



Any linear system with two input and two output chan- 
nels can be described in terms of a 2 x 2 transfer ma- 
trix ^ ^. In many instances (e.g., in polarization op- 
tics H), we are interested in the transformation proper- 
ties of quotients of variables rather than on the variables 
themselves. Then, the transfer matrix induces a bilinear 
transformation that represents the action of the system 
in the complex plane Q . 

For layered media the use of 2 x 2 matrix methods is 
standard, but the corresponding complex-plane represen- 
tation is not a common tool. In fact, it has been recently 
established that for a lossless multilayer the transfer ma- 
trix is an element of the group SU(1,1) which is the 
key for a deeper understanding of its behavior and simpli- 
fies exceedingly its description in terms of bilinear trans- 
formations. 

Moreover, as we have pointed out recently the Iwa- 
sawa decomposition provides a remarkable factorization 
of the transfer matrix representing any multilayer (no 
matter how complicated it could be) as the product of 
three matrices of simple interpretation. At the geomet- 
rical level, such a decomposition translates directly into 
three actions that are the basic bricks from which any 
multilayer action is built. 

In this Letter we go one step further and show that the 
trace of the transfer matrix allows for a classification of 
multilayers in three disjoint classes with properties very 
close to those appearing in the Iwasawa decomposition. 

To maintain the discussion as self contained as possible 
we briefly summarize the essential ingredients of multi- 
layer optics we shall need for our purposes. The configu- 
ration we consider is a stratified structure that consists of 
a stack of plane-parallel layers sandwiched between two 
semi-infinite ambient (a) and substrate (s) media, which 
we shall assume to be identical, since this is the common 
experimental case. Hereafter all the media are supposed 
to be lossless, homogeneous, and isotropic. 

We assume an incident monochromatic, linearly polar- 
ized plane wave from the ambient, which makes an angle 
9o with the normal to the first interface and has ampli- 



tude Ei^K The electric field is either in the plane of 
incidence {p polarization) or perpendicular to the plane 
of incidence (s polarization). We consider as well another 
plane wave of the same frequency and polarization, and 
with amplitude Eg \ incident from the substrate at the 
same angle Oq 0. 

As a result of multiple reflections in all the interfaces, 
we have a backward-traveling plane wave in the ambient, 
denoted E^ •* , and a forward-traveling plane wave in the 



substrate, denoted E^^K 
amplitudes as the vector 



E 



It is useful to treat the field 



Ei+) 

E(-) 



(1) 



which applies to both ambient and substrate media. 
Then, the amplitudes of the fields at each side of the 
multilayer are related by a 2 x 2 complex matrix Mas, we 
shall call a multilayer transfer matrix H , in the form 



where, in our case, 



las can be expressed as 



1/Tas rijt:, 

Ras/Tas 1/^as 
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with det Mas = +1, which shows that the set of lossless 
multilayer matrices reduces to the group SU(1,1), whose 
elements depend on three independent real parameters. 
In this matrix Ras and Tas are, respectively, the overall 
refiection and transmission coefficients for a wave inci- 
dent from the ambient. 

For the case at hand, the Iwasawa decomposition reads 
as i 



where 



Mas = K(,^)A(e)N(i.), 



exp(i0/2) 

exp(-i(/)/2) 



(4) 
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A(0 



cosh(^/2) isinh(^/2) 
-isinh(^/2) cosh(^/2) 

1-ZI//2 v/2 
v/2 l + w/2 



(5) 



In other words, K((/i) belongs to a compact, A(^) to 
an abelian, and U{v) to a nilpotent subgroup, respec- 
tively, and their physical meaning has been discussed in 
Ref. . Furthermore, such a decomposition is global and 
unique 

As we have said, we are often interested in the trans- 
formation properties of field quotients rather than the 
fields themselves. Therefore, it seems natural to consider 
the complex numbers 



(6) 



for both ambient and substrate. From a geometrical 
viewpoint, Eq. defines a transformation of the com- 
plex plane C, mapping the point into the point Zq, 
according to 



a + (3zs 



(7) 



Thus, the action of the multilayer can be seen as a func- 
tion Za = f{zs) that can be appropriately called the mul- 
tilayer transfer function ^ . The action of the inverse 
matrix is Zg = <i>[M^g^, Za]- 

We can define an action of the group SU(1,1) on the 
complex plane C by these bilinear transformations. The 
complex plane appears then decomposed in three regions 
that remain invariant under the action of the group: the 
unit disc, its boundary and the external region ||l0|| . 

The Iwasawa decomposition has an immediate transla- 
tion in this geometrical framework, and one is led to treat 
separately the action of each one of the matrices appear- 
ing in this decomposition. To this end, it is worth noting 
that the group SU(1,1) we are considering appears al- 
ways as a group of transformations of the complex plane. 
The concept of orbit is especially appropriate for obtain- 
ing an intuitive picture of the corresponding action. We 
recall that, given a point z, its orbit is the set of points z' 
obtained from z by the action of all the elements of the 
group. In Fig. 1 we have plotted some orbits for each one 
of the subgroups of matrices K(</)), A(^), and ^{v). For 
matrices K((/)) the orbits are circumferences centered at 
the origin, for A(^), they are arcs of circumference going 
from the point +i to the point —i through z. Finally, for 
the matrices N(i') the orbits are circumferences passing 
through the point +i and joining the points z and —z*. 

To go beyond this geometrical picture of multilay- 
ers, let us introduce the following classification: a ma- 
trix is of class K when [Tr(Mas)]^ < 4, is of class A 
when [Tr(Mas)]^ > 4, and finally is of class N when 
[Tr(Mas)]^ = 4. To gain insight into this classification. 



let us also introduce the fixed points of a transfer ma- 
trix as the points in the complex plane that are invariant 
under the action of Mas; i-e., 



whose solutions are given by 



ilm a ± ^ (Re aY — 1 



(8) 



(9) 



Since Tr(Mas) = 2Re a, one easily check that the ma- 
trices of class K have two fixed points, one inside and 
other outside the unit disc, both related by an inversion; 
the matrices of class A have two fixed points both on the 
boundary of the unit disc and, finally, the matrices class 
TV have only one (double) fixed point on the boundary of 
the unit disc. 

Now the origin of the notation for the classes should be 
clear: if one consider the Iwasawa decomposition (^, one 
can see that the matrices K(<^) are of the class K with 
the origin as the fixed point in the unit disc, matrices 
A(^) are of the class A with fixed points and —i and 
matrices N(^) are of the class N with the double fixed 
point Of course, this is in agreement with the orbits 
in Fig. 1. 

To proceed further let us note that by taking the con- 
jugate of with any matrix C G SU(1,1) we obtain 
another multilayer matrix; i.e.. 



C 



(10) 



such that Tr(Mas) = Tr(Mas). The fixed points of Mas 
are then the image by C of the fixed points of Mas- For 
our classification viewpoint it is essential to remark that 
if a multilayer has a transfer matrix in the class A, 
or iV, one can always find a family of matrices C such 
that the conjugate through C is of the form K(0), A(^) 
or N(j/), respectively. The explicit construction of this 
family of matrices is easy: it suffices to impose that C 
transforms the fixed points of Mas into the corresponding 
fixed points of K(0), A(^), or N(j/). For example, if Mas 
is in the class K and its fixed point inside the unit disc 
is z/, one should have 

$[CMasC-\ 0] = $[CMas, Zi\ - $[C, Z/] = 0. (11) 

In order to appreciate the physical meaning of this clas- 
sification let us take another conjugate of Eq. ([l^) by the 
unitary matrix 



1 

\/2 



1 i 
i 1 



(12) 



In consequence, we can rewrite it alternatively as 

£a=A^as^s, (13) 
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where the new field vectors are defined as 

Ei+) 



£{+) 




(14) 



and the transformed multilayer matrix reads as 

Mas=U MasU-' =UCMasC-^ll-\ (15) 



One can easily check that detAl as = +1 and all its ele- 
ments are real numbers. Therefore, Mas belongs to the 
group SL(2,IR) that underlies the structure of the cele- 
brated ABCD law in first-order optics [|| ^ @. 

By transforming by U the Iwasawa decomposition (^) , 
we get the corresponding one for SL(2,R), which has been 
previously worked out |Tq] : 



Mas = IC{(j>)A{m{i^) 



(16) 



where 



/C(0) 



cos((/)/2) sin(0/2) 
-sin(0/2) cos(0/2) 

exp(^/2) 
exp(-e/2) 

1 0" 

V 1 



(17) 



The physical action of these matrices is clear. Let us 
consider by the moment all of them as ABCD matrices 
in geometrical optics that apply to position x and mo- 
mentum p (direction) coordinates of a ray in a transverse 
plane. These are the natural phase-space variables of ray 
optics. Then K,{(j)) would represent a rotation in these 
variables, -4(^) a magnifier that scales x up to the factor 
771 — exp(^/2) and p down by the same factor, and N{v) 
the action of a lens of power u [Q . 

In the multilayer picture, f^"*"^ can be seen as the corre- 
sponding x, while £^^^ can be seen as the corresponding 
p. Then, the key result of this Letter is that when the 
multilayer transfer matrix has [Tr(Mas)]^ lesser, greater 
or equal to 4 one can find in a direct way a family of 
matrices that gives a new vector basis such that the ac- 
tion of the multilayer, when viewed in such a basis, is 
exclusively rotationlike, or magnifierlike, or lenslike. 

We expect that the formalism presented here could 
provide a general tool for analyzing and classifying the 
multilayer performance in an elegant and concise way. 
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FIG. 1: Plot of several orbits in the unit disc of the elements 
of the Iwasawa decomposition K(^), A(^), and N(i/) for the 
group of multilayer transfer matrices. 



